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Abstract 

We present a summarizing account of a series of inves- 
tigations whose central topic is to address the ques- 
tion whether atomic stabilization exists in an ana- 
lytical way. We provide new aspects on several is- 
sues of the matter in the theoretical context when 
the dynamics is described by the Stark Hamiltonian. 
The main outcome of these studies is that the gov- 
erning parameters for this phenomenon are the total 
classical momentum transfer and the total classical 
displacement. Whenever these two quantities vanish, 
asymptotically weak stabilization does exist. For all 
other situations we did not find any evidence for sta- 
bilization. We found no evidence that strong stabi- 
lization might occur. Our results agree qualitatively 
with the existing experimental findings. 

1 Introduction 

Due to the breakdown of standard perturbation the- 
ory, the understanding of the physics of an atom in 
a strong (intensities larger than 3.5 x 10 16 Wcm~ 2 
for typical frequencies) laser field is still poorly un- 
derstood to a very large extent. Hitherto the large 
majority of the obtained results is based on numerical 
treatments. In our investigations we aim at a rigor- 
ous analytical description of phenomena occurring in 
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this regime. This proceeding will provide an account 
of a series of publications [1-4] . In many cases we will 
simply summarize and state some of the results and 
refer the reader for a detailed derivation to the origi- 
nal manuscripts, but we shall try to put an emphasis 
on new aspects for which we will supply an extensive 
discussion. Several of the presented arguments and 
results may not be found in [1-4]. 

Amongst the phenomena occurring in the high in- 
tensity regime in particular the one of so-called sta- 
bilization has recently caused some controversy, not 
only concerning its definition, but even its very exis- 
tence altogether [5-26]. Roughly speaking stabiliza- 
tion means that atomic bound states become resistent 
to ionization in ultra-intense laser fields. A more pre- 
cise definition may be found in section 2.2. 

2 Physical Framework 

The object of our investigations is an atom in the 
presence of a sufficiently intense^ laser field, which 
may be described in the non-relativistic regime by the 
time-dependent Schrodinger equation in the dipole 
approximation 
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We use atomic units throughout. The time depen- 
dent external electric field will be treated classically 

1 Sufficiently refers to the validity of a classical treatment 
of the laser field. A rigerous quantum electrodynamical treat- 
ment of ionization phenomena was recently initiated in [ p?| . 
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and is assumed to be linearly polarized of the gen- 
eral form E(t) = E f(t), where E denotes the field 
amplitude and f(t) is some arbitrary function which 
equals zero for t < and t > r, such that r defines 
the pulse length. Depending on the context it is con- 
venient to express the Hamiltonian in equation ([!]) in 
different gauges. 

2.1 Gauge equivalent Hamiltonians 

Taking A,<_j (t) to be a one parameter family of uni- 
tary operators, we may construct the gauge equiv- 
alent Hamiltonian Hi (t) from Hj (t) by the usual 
gauge transformation 

H t (t) = idtAj^Ajltit) + A^H^AjUt). 

(2) 

Choosing the most conventional gauge, the so-called 
length gauge, the Hamiltonian to describe the above 
mentioned physical situation is the Stark Hamilto- 
nian 

fff(t) = iff + V(x)+z •£(*). (3) 

V (x) is the atomic potential and = p 2 /2 denotes 
the Hamilton operator of the free particle. We in- 
troduced here sub-and superscripts in order to keep 
track of the particular gauge we are in and to identify 
a specific Hamiltonian, respectively. In our conven- 
tions p and x denote operators, whilst p and x are 
elements in R 3 . Other commonly used Hamiltonians 
are the one in the velocity gauge 

H^t) = \{p-b{t)e z f + V{i) (4) 

and the one in the Kramers-Henneberger (KH) frame 

n 

H$c H (t) = H? + V(5c-c(t)e z ). (5) 

Here e z denotes the unit vector in the z-direction. 
These Hamiltonians may be obtained from each other 
by using 

A v ^(t) = e lb « z , (6) 
A v ^ KH {t) = e - M ('V c «P*, (7) 

Al^ KH {t) = e -ia,{t) e -ib{t)* e ic{t) V% (g) 



in (g). p z is the component of the momentum opera- 
tor in the z-direction. We have employed the impor- 
tant quantities 

t 

b(t) = JdsE{s), (9) 
o 

t 

c(t) = J dsb{s) , (10) 



i 

o(t) = \f dsb 2 (s), (11) 



which are the classical momentum transfer, the clas- 
sical displacement and the classical energy transfer, 
respectively. It will turn out, that in particular b (t) 
and c (r) are the crucial parameters for the descrip- 
tion of the phenomenon we are going to discuss. The 
classical energy transfer a(t) is not a crucial quantity 
since it enters all expressions only as a phase and will 
therefore cancel in all relevant physical expressions. 

In our considerations we will also need the Hamil- 
tonians 

H l A (t) = H?+V(x), Hf v (t) = H?+z-E(t) (12) 

which describe an electron in the atomic potential 
or in the electric field, respectively. Of course these 
Hamiltonians may also be transformed into the other 
gauges by ©-(§). Notice that H% v H {t) = H°. 

2.2 Definition of Stabilization 

Since stabilization means different things to different 
authors and a universally accepted concept does not 
seem to exist yet, we will precisely state our defini- 
tion of it. We will not discuss the behaviour of ion- 
ization rates as some authors do, but we shall con- 
sider exclusively ionization probabilities. Denoting 
by \\ip\\ 2 = = / \ip(x)\ 2 d 3 x the usual Hilbert 

space norm, the ionization probability is defined as 

Vty) = 1-||P + ^|| 2 • (13) 

We used the scattering matrix 

S = lim exp(it + H + ) ■ U(t +: tJ) ■ exp(— it_H_) , 

t± — >±oo 

(14) 
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where H± = lm\ t ^± 00 H(t), ip is a normalized bound 
state of P± is the projector onto the bound state 
space of H± and U(t,t') is the time evolution op- 
erator from time t' to t associated^ to Hit). The 
time evolution operators may be transformed from 
one gauge to another by 

U?(t,t') = A^ t (t)U«(t,t>)Ajl t (t>). (15) 

The ionization probability V(ip) is a gauge invari- 
ant quantity |Q. Note that for the gauge equivalent 
Hamiltonians quoted above we have in general 

lim Hf(t)? lim Hf!(t)? lim H s KH {t). (16) 

r — >±oo t — >±oo r — >±oo 

However, (recall that 6(0) = c (0) = 0), equality in 
the first case holds whenever we have b (r) =0 and 
in both cases when in addition c(t) = 0. We will 
encounter this condition of a particularly switched 
on and off pulse below once more as the necessary 
condition for the presence of what we refer to as 
asymptotically weak stabilization We would like 
to point out that this condition does not coincide 
necessarily with the notion of adiabatically switched 
on and off pulses, because we may achieve b (t) = 
and c (t) =0 of course also with a very rapid switch 
on and off. Since we are interested in the behaviour 
of the atomic bound states \ip(t = 0)) = of the 
Hamiltonian Hf we should commence the discus- 
sion in ([l3]) in the length gauge (in this case we have 
limt^ioo -fff(i) = Hj 4 ) such that in our situation 
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P- 



■Uf (r,0)V| 



(17) 



Regarding the ionization probability as a function of 
the field amplitude Eq, stabilization means that 



dE 



< 



for V (V>) + 1 



(18) 



for Eg G [0, 00) on a finite interval. Hence the oc- 
curence of a saddle point does not qualify as stabi- 
lization. Also we would like to introduce some ter- 
minology in order to distinguish in dl8| ) between the 

2 Associated is to be understood in the sense that the 
time evolution operator obeys the Schrodinger equation 
id t U(t,t) = H(t)U(t,t) . 

3 There are doubts expressed by experimentalists about the 
possibility to realise pulses having simultaneously b (r) = 
and c(r) = & 



case of equality and strict inequality. If the former 
sign holds we call this behaviour "weak stabilization" 
and in the latter case "strong stabilization". In case 
weak stabilization only occurs in the limit E — > oo, 
we shall refer to it as "asymptotically weak stabiliza- 
tion". 



3 Upper and lower bounds for 
the Ionization Probability 

The outcome of every theoretical investigation will 
attach some sort of error to any physical quantity. 
In the minority of cases this error can be precisely 
stated, since it may either be the consequence of var- 
ious qualitative assumptions based on some physical 
reasoning which arc difficult to quantify or it may be 
of a more technical nature originating in the method 
used. For instance for the physical quantity we are in- 
terested in, the ionization probability P(ip), the most 
fundamental error is introduced by the assumptions 
for the validity of the main physical framework, that 
is the Schrodinger equation (|l|) (i.e. non-relativity, 
dipole approximation, classical treatment of the ex- 
ternal field, neglect of the magnetic field, etc.). Ex- 
amples for errors rooted in a particular method used 
are: When solving the Schrodinger equation numeri- 
cally one is forced to discretise H(t), insert the atom 
into a finite box and introduce absorbing mask func- 
tions at the boundary, etc. Also one is not able to 
project on all bound states or all states of the "dis- 
crete continuum" [|] and is forced to introduce a cut- 
off, whose effect is in our opinion not discussed in 
the literature. Some further examples are the errors 
resulting from the termination of a Floquet, Fourier 
expansion or perturbation series^]. We would like the 
reader to keep these basic facts in mind, i.e. "ex- 
act" results do not exist and one is always dealing 
with some form of bounds, when judging about the 
method presented in this section. The essence of the 
method consists in treating bounds which restrict a 
physical quantity rather than looking at its actual 



4 See the conclusion for a discussion of this point. 
5 See section 6 and the conclusion for a discussion of this 
point. 
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value. One of the main virtues of this approach is 
that it may be carried out purely analytically. In 
different contexts it has turned out to be extremely 
fruitful, for instance in the proof of the stability of 
matter pCj and the stability of matter in a magnetic 



field |3J 

We will provide rigorous analytic expressions for 
the upper and lower bound, P u (iJj) and Pi(ip), re- 
spectively, for the ionization probability in the sense 
that 



PiW<P(*P)<PuW 



(19) 



Hence within the basic theoretical framework upper 
and lower bounds serve as sharp error bars. Surely 
one should treat these expressions with care and be 
aware of their limitations in the sense that about the 
actual shape of P(ip) no decisive conclusion can be 
drawn whenever Pi{ip) differs strongly from P u (ip). 
However, it seems a reasonable assumption that the 
analytic expression of the bounds reflect qualitatively 
the behaviour of the precise ionization probability. 
Nonetheless, there exist certain type of questions in 
the present context which can be answered decisively 
with this method. Concerning the question of stabi- 
lization we may consider the bounds as functions of 
the field amplitude and can conclude that stabiliza- 
tion exists or does not exist once we find that Puiip) 
for increasing field amplitude tends to zero and Pi(if)) 
tends to one, respectively. Unfortunately, one does 
not always succeed in deriving analytic expressions 
which are of this restrictive form. 

In Q we obtained 



PM) = 1 - I J \\(V(x-c(t)e z ) - V{x))iA\dt 
o 

2 



2E + 6(r) 2 



2\b(r) 



2E + b(r) 2 



\{V(x-c{T)e z )-V(x))i>\\ 



\pzM\ 



(20) 



which is valid when — E < 6(r) 2 /2. Here E is the 
binding energy. With the same restriction on 6(r) we 



found as an upper bound 



Pu(ip) = < \\{V(x-c{t)e z )-V{x))^\\dt (21) 



+ |c(r)||M.|| + ^W M / ' 

By a slightly different analysis we also derived a 
bound valid without any additional restrictions 



PuW = < J \\(V(x-c(t)e z )-V(x))ip\\dt 
o 

+|c(r)| \\p^\\ + \b(r)\ . (22) 

In |2| we applied these bounds to the Hydrogen atom, 
obtaining 

Wn00) = l- 



j3/2 5( r )2 _ l/ n 2„3/2 

1 2|6(r)| I 2 



(23) 



PuO'nOo) 



2r \c(t)\ 5n 3 



i 3 / 2 nV3 



G 



'■Hr)\ (24) 



Our method allows in principle to consider any bound 
state, but initially we restricted ourselves to s-wave 
functions, keeping however the dependence on the 
principal quantum number n. The discussion of ( |2l| ) 
and ( ^ ) in [jlj was plagued by the requirement that 
the pulse duration should be fairly small. This limi- 
tation, which ensured that bounds for the ionization 
probability are between zero and one, was overcome 
in [^), since in there we had an additional parameter, 
i.e. n, at hand. As the expressions ( p3| ) and ( |2^ ) 
show after a quick inspection, one may achieve that 
their values remain physical, even if one increases the 
pulse duration, but now together with n. 

In particular we investigated the effect resulting 
from different pulse shapes, since it is widely claimed 
in the literature that a necessary condition for the 
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existence of stabilization is an adiabatically smooth 
turn on (sometimes also off) of the laser field. For 
definiteness we assumed the laser light to be of the 
general form E(t) = Eos'm(uit)g(t). Besides other 
pulses we investigated in particular the ones which 
are widely used in the literature, where the envelop- 
ing function is either trapezoidal 



C ± for 
g(t) = { 1 for 
1 ^ for 

or of sine-squared shape 



< t < T 

T < t < (t — T) (25) 
(r-T)<t<T 



!sin 2 (g,) for < t < T 

1 for T<t<(r-T) 

sin 2 (4^) for (r-T)<t<T. 

(26) 

We found no evidence for stabilization for an adi- 
abatically smoothly switched on field, provided that 
b{r) 7^ 0. The latter restriction emerges in our anal- 
ysis as a technical requirement, but it will turn out 
that it is of a deeper physical nature^] 

4 Ionization Probability in the 
ultra-extreme Intensity Limit 

Since we are interested in very high intensities we ex- 
pect to be able to draw some conclusions from the ex- 
pressions for the ionization probability in which the 
field amplitude is taken to its ultra-extreme limit, 
i.e. infinity. In particular we may decide whether 
asymptotic stabilization exists. Despite the fact that 
in this regime one should commence with a relativis- 
tic treatment, our physical framework, that is the 
Schrodinger equation (|l|) remains self-consistent, and 
should certainly represent the overall behaviour. In 
we rigorously take this limit under certain gen- 
eral assumptions on the atomic potential^ and the 

6 Already in (|lrj|) we observed that fe(r) = is somewhat 
special. 

7 We made the general assumptions that V (x) is a real mea- 
surable function on R 3 and that each e > one may decom- 
pose V as V = V\ + V2 where Vi is in L 2 [R 3 ) (i.e. square 



laser field, which include almost all physical situa- 
tions discussed in the literature. Whenever 6(r) and 
c(t) vanish simultaneously we found 

lim V (ip) = 

whereas in all other cases we obtained 

lim V (V>) = 1 . (28) 

|-Eo|^oc 

This means that in the former case we have asymp- 
totically weak stabilization. It should be noted that 
in this analysis the pulse shape is kept fixed, such that 
adiabaticity can not be guaranteed anymore. Hence, 
weak stabilization is found for a situation in which 
it is generally not expected to occur. It would be 
very interesting to perform similar computations as 
in H in which the pulse shape is varied in order to 
keep also adiabaticity and study whether the effect 
will become enhanced in any way. Furthermore, in 
this case the time evolution operator coincides with 
the one of the free particle H® 

lim II (Uf (r, 0) - exp -itH^) ip\\ = 0. (29) 

|E |— >oo 

Of course this type of argument does not allow to 
draw any decisive conclusions concerning strong sta- 
bilization. 

5 Gordon- Volkov Perturbation 
Theory 

In the high intensity regime for the radiation fields, 
the basic assumption for the validity of conventional 
perturbation theory breaks down, i.e. that the abso- 
lute value of the potential is large in comparison with 
the absolute value of the field. However, there is a re- 
placement for this, the so-called Gordon- Volkov (GV) 
perturbation theory |}2| . Since the basic idea is sim- 
ple, it makes this approach very attractive. Instead 

integrablc) with compact support and V2 is in L°° (H 3 ) with 

Halloo =ess SU P |V2(a?)| < £• Furthermore we assumed that 

SG-R 3 

has no positive bound states. Such potentials are Kato 
small. In particular the Coulomb potential is Kato small. 



-irH) 



< 1 



(27) 
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of constructing the power series, either for the fields 
or for the time evolution operator, out of the solution 
for the Schrodinger equation involving the Hamilto- 
nian Hf- and regarding zE(t) as the perturbation, 
one constructs the series out of solutions involving 
the Hamiltonian Hf v and treats the potential V as 
the perturbation. 

The starting point in this analysis is the Du Hamel 
formula, which gives a relation between two time evo- 
lution operators U°"(t,t') and f/,.(t, t') associated to 
two different Hamiltonians Hf(t) and Hj(t), respec- 
tively 

U?(t,t') = U*(t,t')-if ds U? (t,s)Hlf(s)U^s,t'). 
t> 

(30) 

Here we use the notation H^- (s) = Hf(s) — H^(s). 
The formal iteration of (|3^) yields the perturbative 
series 

oo 

^(M') = EO"M') ■ ( 31 ) 

71=0 

We introduced in an obvious notation the quan- 
tity (n\t, t') relating to the time evolution 
operator order by order in perturbation theory, 
i.e. U*j>(0\t,t') = U$(t,t% U«f(l\t,t') = 

i f*ds uf(t,s)Hlf(s)U^(s,t'), etc. It should be 
noted that the perturbative series is gauge invariant 
in each order, since U^j{n\t, t') is a gauge invariant 
quantity by itself. Mixing however expansions for 
different choices of the Hamiltonians, i.e. a and b 
or different gauges i and j will not guarantee this 
property in general. A rather unnatural choice (for 
instance with regard to the possible convergence of 
the series) would be i ^ j. Taking therefore i = j 
and in addition a = S and b = GV we obtain 

Uf(t, t) = ur(t, t) + U*f v {l\t, t') + ... (32) 

In this case we need H?f V ' (t) = H^ v (t) = V (x) , 

H^ H V KH (t) = V(x — c(t)e z ) and the Gordon- Volkov 
time evolution operator, which in the KH-gauge 
equals the free-particle evolution operator in the 



length gauge 

U^(t,t') = A^ KH (t)ur(t,t')A^ KH (t<) 

= U?(t,t'). (33) 

The expressions for the Gordon- Volkov time evolu- 
tion operator in the length and velocity gauge may 
then simply be obtained from d33| ) by the applica- 
tion of (@)-(^) according to (|l5|). The choice i = j 
together with a = S and b — A in (^) yields the 
usual perturbation series, which is well known from 
the low intensity regime. One may also decide for a 
rather strange procedure and take the latter choice 
in the first iterative step and terminate the series af- 
ter the second iterative step in which one makes the 
former choice. In that case one obtains 

Uf(t,t) = U z A (t,t) (34) 
t 

i J ds U A {t 1 s)H?f{s)U? v {s,t') + 0(n 2 ). 
v 

For i = I or i = v, this procedure is sometimes re- 
ferred to as the Keldysh- Q or Faisal-Reiss [|4| ap- 
proximation, respectively. As we demonstrated this 
method is of course not "non-perturbative" , as some- 
times wrongly stated in the literature. 

There are some exact results which may be derived 
from the perturbative expression, one concerning the 
ultra-extreme intensity limit of the previous section 
and the other the ultra-extreme high frequency limit. 
Both results are simple consequences of the Riemann- 
Lebesgue theorem^. We obtain 

lim Ai^At) = 1 (35) 

u — >oo 

such that with j3^ ) 

lim U^(t,t') = U?{t,t') = e-^-^ H °. (36) 

UJ — >oo 

Since the atomic potential is independent of Eq we 
obtain with ([36]) that the entire Gordon- Volkov series 

8 If g(x) E L\ (—00,00) (i.e. |ff(^)| is intcgrablo) then 

00 

lim I g{x)e~ ltx dx = 0. 

t— >±oo J 
— 00 
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( pl| ) is independent of the field amplitude as well, 
such that 



4r ( lim 







(37) 



We have therefore weak stabilization in this ultra- 
extreme high frequency limit for all systems for which 
( pl| ) makes sense and for which the laser field is of the 
form E(t) = Eq sm(ojt)g(t). One should keep in mind 
that the use of the Stark Hamiltonian in ([j]) assumes 
the validity of the dipole approximation, such that 
the limit u — > oo (a slightly milder assumption was 
used in the seminal paper || which formulates the 
high frequency approach) only makes formally sense. 
In order to describe real physics in this frequency 
regime one should actually also take multipole terms 
into account. This is a further example of the errors 
we mentioned at the beginning of section 3. 

Concerning the ultra-extreme intensity limit, we 
consider the transition amplitude between two bound 
states ipt (x) , tpj (x) of the Hamiltonian Hf perturba- 
tively 



<^,C/f(r,0)V^) 

= (i; l ,A KH ^ l (T)U^ H (T,0)^ J ) 

= {l> i ,A KH ^t(T)U$%(T,0)ll> j ) 

+ (A,Akh^i(t)U%%V kh {1\t,0)1> 3 



(38) 



Recall that a(0) = b{0) = 
Ai^kh{0) = 1. Using now 
zeroth order we obtain 



c(0) = 0, such that 
fc33h it is clear that to 



,-iVfl) 



tpJx) 



(39) 



when b(r) = c(r) = 0. In all other cases we may bring 
this term into a form suitable for the application of 
the Riemann-Lebesgue theorem, such that the zeroth 
order matrix element always vanishes in the ultra- 
extreme intensity limit. For the higher order terms 
the argument is analogous with the difference that 
the condition b(r) = c(r) = does not have the con- 
sequence that these expressions become independent 
of Eq, since also terms like b(t),c(t) for < t < t 
appear. Hence by the application of the Riemann- 
Lebesgue theorem all higher order terms vanish in 



the limit Eq — > oo. If we now sum over all bound 
states i in ( |38| ) we obtain the results of section 4 ( |27| ) 
and (p8|). 

6 1-dimensional potential 

In Q we applied the GV-perturbation theory to the 
one-dimensional delta potential with coupling con- 
stant a 

V(x) = -aS(x) . (40) 

In the momentum space representation this potential 
becomes 

V(p,p>) = (p\V\pi)=~^- (41) 

and the wave function for its only bound-state is well 
known to be 



a- 2 



7r a 2 + p 2 



(42) 



21 



The very fact that this potential possesses one bound 
state only with bound state energy —a 2 /2 makes 
it a very attractive theoretical atomic toy potential 
(e.g. 



, f36|). Also with regard to the GV- 
perturbation theory one expects intuitively a good 
convergence. 

We construct the exact time-dependent wave func- 
tion: 

ip (p, t) = 4> GV (p, t) + * ( P , t) 

with 

t 

*(p,t) = i— f dse-^We^M*)) 
2tt J 



(43) 



xe 2 



(p-Ht)f(t-s) 



1>i is) 



4>i H) 



W / dpil>(p,t) 



(44) 
(45) 



Integrating (43) with respect to p we obtain a 
Volterra equation of the second kind in t 



it) 



ice^ I e 2 <*~ s ) 

dpip G v(p>t)+\l J dsipj (s)— ==. 



(46) 
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The virtue of this equation is that the error of its 
solution, even when obtained by an iterative proce- 
dure, is completely controllable. The iteration of the 
Volterra equation yields 

oo 

« oo 
t/>(t)= J dp^ GV {p,t)+Y,^n{t) (47) 



with ip n (t) denoting the function order by order. We 
derived Q an upper bound for the absolute value of 
this function 

\i> n (t)\ < V8^-^ (aVt/2) n • (48) 

It is known from the theory of integral equations fl3~7| 
that this is sufficient to prove that the series con- 
verges for all values of a, and in addition we were 
even able to sum up the whole series 



| = V^(2exp(a 2 i/2)-l 

n=l 

-U U {a 2 t/2) /V?F) (49) 

such that we can compute the maximal relative error 
after the zeroth order in GV perturbation theory 



fl = 2v / 7T 



(2exp(a 2 t/2) - 1 - Ux^ (aH/2) /0F 



(Ul 1 (*_) +E7l 1 ($+) 

\ 2 > 2 x ' 2 > 2 v ' 



Here Ui i (z) is the confluent hypergeometric func- 
tion (see for instance |3q]) and 



$± :=tc? (±7+-(l-7 2 ) ) ,7 



c(r) 



ra 



(51) 



This analysis allows us to determine the error which 
is introduced by the termination of the GV-series. 
Usually this is done without any justification about 
the precision and the only reasoning provided is very 
often solely the comparison with the next order term. 
This is however not enough as one knows from simple 
rest term estimations in a Taylor expansion. 
The ionization probability turns out to be 



W)Y 



(52) 



= l-|(V;,V G yM) + (^(T)>| 2 (53) 



with 



(^,V G y(r))--a 3 e- M M 



(54) 



dp 



exp ( — ir^; — ic (r) p 



(a 2 + (p + 6(r)) 2 ) {a 2 +p 2 ) 



and 



(V,*(t)) = i 



a 5 f 

—tjds^jis) 





(55) 



dp- 



a ic TS (p-b(r)) „-f(p-b(-r)) (t-s) 



(a 2 +p 2 ) 



\q(T)\ has of course the interpretation as survival 
probability. We use the abreviation c U ' '■= c(t) — 
c(t'). These expressions constitute explicit examples 
for the general statements made at the end of section 
5. Taking the ultra-extreme intensity limit we obtain 
as a consequence of the Riemann-Lcbesgue theorem 



lim \q(r)\' 

Eq— »oc 



\{^av(rW for b(r)=c(T) 
otherwise 







(56) 

This means we have asymptotically weak stabiliza- 
tion and the result is in agreement with the one in sec- 
tion 4. However, the assumption made on the poten- 
tial in H does not include the 5-potential[], such that 
( |5^ ) is not only obtained by an alternative method 
but also covers an additional case. As we already 
observed, a different physical behaviour is obtained 
depending on the values of b(r) and c(r) and it is 
therefore instructive to treat several cases separately: 
To the lowest order we obtained 

i) 6(r) = 0, c(t) = 



V (V) = 1 - - 

7T 



Ui 



(57) 



3 The <5-potential is not a Kato small potential. 
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ii) 6(r) = 0, c(t) ^ 

g(r) = p_ Uii ($+) + I7i i ($_) 



-ira {/ 1 3 



(58) 



with yj ± = | (l ± ac(r) — ira 2 ). Equation 
is only valid for |7| < 1. 
iii) b(r) yt 0, c(r) 



?(t) = ^| C/i i($+)+C/i 



(59) 



P U ($-) + [/u($ + ) 



with 



6(t) 2 ± 2a6(r)i 



$± = ra 2 ( ±7 + - (l - 7 2 ) 



(60) 
(61) 



Equation ( |59| ) is only valid for |7| = 



< 



1. 

The restrictions on the parameters 7 and 7 orig- 
inate in the limited validity of the integral repre- 
sentation for the confluent hypergeometric functions, 
which is employed here. To obtain all higher orders 
we have to compute 



(V,*(t)) = 



(62) 



2ttV2 

ds^jj (s) e^T (f/x 1 ($'_) + ($' + )) 

±a(ct s -(t-s)b(t)) 

+|(*-) ^-(^- 6 (t)) 2 ).(63) 



Hence the problem to compute the ionization prob- 
ability has been reduced to solving (^6j) and subse- 
quently evaluate (b2h. Surely this is not possible to 



with 



perform entirely in an analytical way, but the initial 
problem has now been reduced to a numerical task, 
whose error is well under control. In Q we carried 
out this analysis for a pulse involving the trapezoidal 
enveloping function. We do not find any evidence for 
strong stabilization even for b(r) = and c(r) = 0, 
however, asymptotically weak stabilization exists for 
the latter case. 



7 Conclusions 

The main outcome of our investigations is that the 
governing parameters for the behaviour of an atom in 
an intense laser field are the total classical momen- 
tum transfer b{r) and the total classical displacement 
c(r)Pl Whenever both these two quantities vanish, 
asymptotically weak stabilization does exist. Also 
the authors of p(| found asymptotically weak stabi- 
lization for this situation (see Fig. 3 therein). For all 
other cases we did not find any evidence for stabiliza- 
tion. 

Since our findings apparently differ from many the- 
oretical results other authors obtained by alternative 
methods we would like comment on possible resolu- 
tions for this discrepancy: 

Introducing a cut-off in the number of bound states 
will produce an upper bound for the ionization prob- 
ability, since in this way one effectively enlarges the 
continuum. In case one finds stabilization for such a 
bound one could confidently conclude that this effect 
indeed exists. Since for lower intensities one can cer- 
tainly expect this bound to be relatively close to the 
real value, whereas for higher intensities this bound 
should decrease even more. 

The introduction of a cut-off (e.g. equation (3) in 
|l2]] ) in the "discrete continuum" (besides the fact 
that this is an ill-defined concept) yields a lower 
bound for the real ionization probability, which is ex- 
pected to be relatively accurate for low intensities 
but very far from the real value for high intensities. 



10 Some authors claim [ f^t that for a pulse of the form ( |25[ ) 
b(T) and c(T) should be relevant parameters. Re sides the fact 
that for a smoothly differentiable pulse (e.g. few)) these quan- 
tities are not precisely defined they do not emerge in our anal- 
ysis as significant. 
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Hence keeping the cut-off fixed and interpreting the 
result obtained in this way as "exact" , one has cer- 
tainly introduced an artificial mechanism to "create" 
stabilization. 

A further approach is based on the Fourier expan- 
sion of the Hamiltonian in the KH-gauge and there- 
after simply keeping the lowest terms. The way this, 
in principle legitimate, method is carried out some- 
times makes conceptual assumptions which are in 
clear conflict to our main physical framework. For 
instance the basic Hamiltonian used in |p2|| (equation 
(1) therein) is in our notation H^ H (t) for an instan- 
taneously switched on monochromatic laser pulse, i.e. 
E(t) = Eq cos(uit). However, the authors of (22) claim 
stabilization to exist for pulses with a smooth adia- 
batic turn on and off (gaufiian enveloping function). 
Clearly equation (6) in p2| breaks the gauge invari- 
ance discussed in section 2.1. of our manuscript, 
such that the authors consider an entirely different 
system. In other words the Hamiltonian (6) in |^2| 
is not gauge equivalent to Hf(t) for a laser pulse 
with gauBian enveloping function. The potential in 
H^- H (t) is shifted by c(t) and not E(t). From our 
point of view the findings of the authors are not sur- 
prising since they artificially impose that c(t) = 0, 
(we also find asymptotically weak stabilization in this 
case) such that one should solve the inverse problem 
in this case to find out which pulse is really consid- 
ered there. However, even under these assumptions 
we would still not find strong stabilization. The same 
procedure is used for instance in (23[ 53. 

Concerning the investigations which do not find 
any evidence for stabilization at all, we would like to 
make the following comment on the use of the GV- 
perturbation theory. In the last reference of @ a 
pulse with instantaneous switch on was used in this 
context, i.e. E(t) — Eq cos(ujt), and an analysis up to 
first order GV-perturbation theory was carried out. 
Typical parameters in were a — 1/2, E = 5, w = 
1.5 and the pulse length was 2 cycles, that is r ~ 8. 
For these parameters we obtain for the relative max- 
imal error /i « 8.44, such that we do not expect the 
GV-perturbation series to be a good approximation 
up to this order and statements made in this context 
should be treated with extreme care. 

We would like to conclude with a remark on the 



existing experimental findings E3j. So far the ex- 
periments carried out only find evidence for asymp- 
totically weak stabilization and confirmations for the 
existence of strong stabilization do not exist to our 
knowledge. We are therefore in complete qualitative 
(it is difficult to determine which values 6(r) and c(r) 
have for the experimentally employed pulses (see also 
footnote 3 for this)) agreement with the existing ex- 
periments. 
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